On the Short Distance Part of the QCD 
Anomaly Contribution to the b — ► sr]' 

Amplitude 



Jan 0. Eeg a '*, Kresimir Kumericki b , Ivica Picek b 

^Department of Physics, University of Oslo, N-0316 Oslo, Norway 

h Department of Physics, Faculty of Science, University of Zagreb, P.O.B. 331, 

HR-10002 Zagreb, Croatia 



Abstract 

In addressing the B — ► rj'K puzzle, there has been a considerable hope in the 
literature to resolve it by the QCD anomaly contribution to the b — > sr]' amplitude. 
This contribution corresponds to the electroweak b — > sg*g* transition followed by 
the off-shell gluon fusion g*g* — ► rj' . In the present paper we perform a critical 
reassessment of this issue. We show that for the hard virtual gluons in a loop there 
is a well defined short distance amplitude corresponding to a remnant of the QCD 
anomaly. However, we find that it cannot account for the measured amplitude. 

In addition, we point out that the reduction of the gluon fusion vertex for the 
off-shell gluons is compensated by an absence of the claimed suppression in the 
electroweak vertex, and that some nonperturbative contributions related to the 
QCD anomaly may still be viable in explaining the physical B — > rj'K amplitude. 
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1 Introduction 



Recent measurements [1,2,3,4,5] of two-body charmless hadronic B meson 
decays confirm the so-called B — > rj'K puzzle. According to the PDG average 
[6] of CLEO, BaBar and Belle measurements, the B — > Krj' decay rate turned 

* Corresponding author. 

Email addresses: j .o.eeg@fys.uio.no (Jan O. Eeg), kkumer@phy.hr (Kresimir 
Kumericki), picek@phy.hr (Ivica Picek). 



Preprint submitted to Elsevier Science 



1 February 2008 



out to be unexpectedly large when compared to the rate obtained within the 
standard effective Hamiltonian approach using quark operators multiplied by 
Wilson coefficients [7]. Namely, the enhancement of 

Br(B + -> K + r]') = (7.5 ± 0.7) • 10" 5 , 
Br (B° -> K°rf) = (5.8+};J) • 1(T 5 , (1) 

when compared to the QCD-penguin dominated B — > Kir rates 

B(B + -> ^°7T + ) = (1.73+J;^) • 10" 5 , 

B(B° -> AT+tt - ) = (1.74 ±0.15) • 10~ 5 . (2) 

calls for an additional contribution to the i^r/' channel, either by performing 
a more complete treatment within the standard model (SM) or by invoking 
new physics beyond the standard model (BSM). Since the properties of the 
//-particle are related to the QCD axial anomaly, it has been quite generally 
expected that the enhancement in (1) is related to this anomaly. 

Among various mechanisms considered to explain the B — > Krj' amplitude 
within the SM, the b — > srj' transition plays a distinguished role. There was an 
invitation in [8] to consider it as a promising short- distance (SD) mechanism. 
It pertains to the recently studied singlet penguin mechanism [9,10], involving 
two gluons fusing into r]'. 

In the present paper we focus on the SD contribution to b — > si]' incorporat- 
ing the electroweak b — > sg*g* transition followed by the off-shell gluon fusion 
9*9* — ¥ rf, where g* denotes a virtual gluon, and rf ~ i] corresponds to the 
flavour singlet state. For low-energy gluons the g*g*f]' vertex should be dom- 
inated by the QCD triangle anomaly analogously to the tt° — > 27 amplitude. 
However, for (at least) one highly virtual gluon of momentum q, the g*g*i] / 
vertex is expected to be suppressed like 1/q 2 similarly to the effective 7*771-° 
vertex [11,12,13]. This remnant of the anomaly for hard off-shell gluons we 
will call the "anomaly tail" contribution. There has been a renewed interest 
in such contributions in view of the hadronic contributions to (g — 2) that they 
induce [14,15]. 

A more complete discussion on the various mechanisms presented in the lit- 
erature for B — > rj'K [16,8,17], including various versions of form factors for 
g*g*T)' [18,19,20,9,10,21,22] will be relegated to a forthcoming paper. In these 
proposals, as a rule, the gluon virtuality employed in the g*g*r)' transition form 
factor lies below mf. In contrast, we consider highly virtual gluons, above the 
nib scale, corresponding to the mentioned anomaly tail contribution. 

After presenting in Sect. 2 the electroweak b — > sg*g* vertex, in Sect. 3 we 
demonstrate how the virtual gluons are glued to the "anomaly tail" part of 
the g*g*rj vertex mentioned above. The resulting contribution is dubbed the 
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Fig. 1. The one-loop graphs which contribute to the box (a) and triangle (b) parts 
of the amplitude (3). 

"short distance anomaly" (SDA) in what follows. In the concluding section 
we discuss the meaning of our results and the relation of our contributions to 
those existing in the literature. 



2 The Flavour Changing b — ► sg*g* Transition 



The flavour changing transitions into two virtual gluons were considered by 
two of the authors in the context of the double-penguin contributions to the 
K° — K° [23] and B° — B° mixing [24]. They have been subsequently studied 
by Simma and Wyler [25] in the case of rare S-decays. Let us now reconsider 
these transitions by taking the symmetric gluon momenta, so that we will be 
able to present the analytical expressions suitable for evaluation of hard-gluon 
loop integrals. The resulting b — > sg*g* amplitude reads 

M^(b^ S g*(p)g*(-p)) = i^^st a 't a £ KH^b + (crossed), (3) 

V 1 71 V 1 i=u,c,t 

where t a denote colour matrices, and Aj are the Cabibbo-Kobayashi-Maskawa 
(CKM) factors. H l subsums the contribution from the box on Fig. 1(a), the 
contribution from the triangle on Fig. 1(b), and the contribution from the off- 
diagonal b — > s self-energy. The b — > s g*(p) loop in Fig. 1(b) is proportional to 
the gluonic monopole, (p 2 7 M — p ■ 7 p M ) L . This p 2 dependence is canceled by 
the 1/p 2 dependence of the s-quark propagator in this reducible diagram. After 
combining the triangle and the self-energy contribution, the UV divergences 
mutually cancel, and only the monopole, triangle part contributes significantly. 

By anticipating the antisymmetric structure of the r]'g*g* vertex into which 
the virtual gluons on Fig. 1 proceed, we select the relevant antisymmetric 
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contribution: 



iff = (—ie^^p^pL) Ai + (/i-v symmetric part) , 
where we obtain 
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The abbreviations are 
D 



Y 1 = l- 



I D 1 
Y2_ D^C ln c"C' 



D = xm\ + (1 — x)M w — x(l — x)p 2 , C = m\ — x(l — x)p 2 . 

This amplitude agrees with that of [25] in the region of their mutual validity. 
Let us stress that Ai has the asymptotic behaviour 1/p 2 for high gluon mo- 
menta (— p 2 — > oo), which is essential in order to obtain an overall finite gluon 
loop contribution to b — > si]' on Fig. 2. Another interesting limit is the leading 
logarithmic approximation 
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where we have neglected the w-quark contribution which is CKM suppressed. 
It should be noted that this leading contribution comes from the triangle 
graph in Fig. 1(b), corresponding to the first term in (5). The box graph has 
no leading logarithm and is numerically small. This clearly differs from the 
dominance of the box part displayed on Fig. 1 of [10], and an explanation on 
the relevance of the b — > sg*g* amplitude (3) is in order. 

There were suggestions in the literature (e. g. [17]) that because of large cancel- 
lations observed by [26,27,25], the contribution of the b — > sg*g* mode followed 
by a g*g* — > rj' transition would be extremely small. It should be strongly em- 
phasized that this statement, valid for soft gluons, is wrong for hard virtual 
gluons, as explicated in the expressions above. Therefore, in the next section 
we take under scrutiny the amplitude stemming from the antisymmetric part 
of the electroweak amplitude on Fig. 1. 
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Fig. 2. The hard gluon loop contribution to the b — > srf transition, determined by 
vertices (3) and (7) 

3 Short Distance Anomaly Contribution 



In this section, we will consider the anomaly tail part of the effective g*g* — > rj' 
vertex to be connected to the b — > sg*g* amplitude from the preceding section. 
By restricting to the hard gluons, having a virtuality above the m& scale, we 
aim at singling out the short distance amplitude for the b — > srf transition. 
This amplitude is of the same order in a s as recently studied amplitudes with 
(colour) singlet penguin topology [28,29,9,10]. 

We need an expression for the gluonic triangle contribution for g*g* — > rf on 
Fig. 3. Using the kinematical choice q\ = K/2 + p, q 2 = K/2 — p for the two 
gluon momenta in this figure, we get a general expression for the g*g* — > rf 
vertex 

N$ /(ft) - rf(K)) = -iS aa 'e^p a K? G(p 2 ). (7) 

In our case, G(p 2 ) will turn out to be a remnant of the gluonic anomaly. A 
priori, the quantity G depends also on the momentum K, but for our purposes 
we only need to keep K to first order there, which means that G is taken as 
independent of the r/ momentum. 

Our effective vertices (3) and (7) allow us a perturbative evaluation of the 
amplitude displayed on Fig. 2: 



AsMb - si) = J Mp' (b s g* g*) (g* g* - r, 



p2 



Having in mind that both vertices M and imply hard gluons in the loop, 
we arrive at the short distance representation of the singlet penguin diagram 
discussed in [28,29]. After substitution of (3) and (7), our result for the short 
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distance di-gluon mechanism reads 



AsMb - si) = 2i-^{sK- 1 Lb) £ \ t J 



d 4 p 1 a s 



MP 2 ) G(p 2 ) . (9) 



(2tt) 4 p 2 7T 



Since Aj(p 2 ) and GQo 2 ) are already one-loop quantities, As D a(6 — > srj') is a 
3-loop amplitude. 

Let us stress that eqs. (8) and (9) apply quite generally for all g* g* rj' form 
factors G(p 2 ). In this paper we will use these formulae for the SDA contribution 
only. 

The evaluation of the g*g* — > rf transition amplitude for off-shell gluons is 
theoretically analogous to the evaluation of the "off-shell anomaly" for the 
flavour triplet axial current related to tt° — > 7*7*. A naive treatment of this 
photonic case leads to problems with unitarity for the e + e~ — > 7* — > 7r°7 
amplitude presented by Jacob and Wu [30], and a number of papers [11,12,13] 
has been devoted to its cure. 

The pertinent triangle amplitude was calculated first by Rosenberg [31] and 
reconsidered by Adler [32], who observed difficulties with divergences when 
inserting the triangle loop into the next loop. On the other hand, for off-shell 
photons, it was shown by [12,13] that the perturbative mass- independent part 
in the triangle diagram is canceled by the pion pole anomalous contribution. 
This observation was essential to resolve the Jacob- Wu paradox [30]. In our 
gluonic case it means that, although the quantity G in (7) obtained in a 
perturbative calculation of the quark-triangle loop is a priori proportional to 



representing the denominator emerging from the triangle-loop integral. 

For the g*g* — > rj transition the flavours j = u,d,s contribute (eventually 
only the s-quark contribution turns out to be important), and we obtain 




(10) 



the unity term gets canceled [12,13]. In (10) we introduced 

Qj = vi 1 - y)ii + x ( x - x )i2 + 2x v qi-12- m] 



(ii) 
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where f n > ~ f n ~ 92 MeV and 
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Fig. 3. The rj — > (7*5* amplitude with indicated symmetric kinematical choice. 



This represents the short distance amplitude for hard gluon virtualities far 
above the r/'-mass. This effective vertex enters into the final gluonic loop on 
Fig. 2 for loop momenta p 2 ^> K 2 = m 2 ,. In the evaluation of the loop integral 
producing the short distance 3-loop contribution for the b — > srj' amplitude, 
the overall renormalization scale \i will be of the order m b . Taking into account 
that /j, acts as the effective IR cut-off in our loop integral, we find that the 
expression (13) is slightly modified (Qj — > Qj + fi 2 in (13) above). In the 
leading logarithmic approximation we obtain for the s-quark term 

'--(3W*)- 

while the u- and ci-quark contributions are completely negligible due to their 
small masses. On account of the CKM unitarity, the loop integral parts in (9) 

= r* ^' <15) 

form the relevant GIM combination displayed on Fig. 4. This plot exhibits 
a reasonably mild dependence on the infra-red cut-off (/x) in the gluon loop, 
justifying a sensible short-distance amplitude. Note that the UV convergence 
of the loop integral is guaranteed by the high energy behaviour of (5) and (14). 
Since the used form-factor (14) corresponds to the anomaly tail term, the net 
contribution can be termed the "short distance anomaly" , as suggested in the 
Introduction. Also, Fig. 4 shows the dominance of the triangle with respect to 
the electroweak box contribution. In the leading logarithmic approximation, 
we obtain from eqs. (6), (14) and (15) 

which may be compared with the result from the (octet) penguin operator to 
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Fig. 4. The infra-red stability check of the anomaly tail contribution (solid line) 
displaying also the dominance of the triangle part (dashed line) over the box part 
(dotted line) from Fig. 1. 



lowest order in the leading logarithmic approximation 

G F 



Mp eng — A c 
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(17) 



4 Discussion and Conclusions 



In the present letter we attempt to clarify a possible role of the b — > si]' 
transition in explaining the B — > rj'K amplitude, and in particular the role 
of the QCD anomaly in obtaining the b — > si]' amplitude. The result above 
shows that we are able to successfully distinguish the short distance b — > si]' 
amplitude (SDA) related to the QCD axial anomaly. Thereby, the involved 
form of the flavour- changing vertex for hard gluons is in clear contrast to 
the result obtained by [10] for soft gluons, where the box (Fig. 1(a)) instead 
of the triangle (Fig. 1(b)) dominates. This enhancement of the electroweak 
vertex for the hard off-shell gluons is compensated by a suppression in the 
two gluon fusion to i]' vertex. The net b — > si]' amplitude coming from this 
anomaly tail turns out to be dominated by the strange quark contribution. In 
the present approach it corresponds to the quark triangle contribution to the 
rj'g*g* coupling for the highly off-shell gluons. We note that the SDA result 
vanishes in the chiral limit m, — > 0. 



8 



Thus, demonstrating that we have obtained a contribution which is different 
from those already existing in the literature, we can attempt to compare it to 
some related amplitudes. A simple comparison of our singlet amplitude (9) at 
the leading logarithm level (16) to the ordinary penguin in (17): 



where we didn't explicate some additional factors in Asda and Ap eng which 
cancel in the ratio. This ratio is at the level of a few percent, but depends 
strongly on what one takes for the involved parameters. How this estimate 
will be modified within a more proper treatment of renormalization group 
equations is also an issue to be further investigated. We have done a simple 
estimate which shows that short distance QCD corrections do not change this 
result substantially. 

To conclude, the above demonstrates that purely short-distance anomalous 
(SDA) aspects of the b — > srj' vertex are marginal in explaining the B — > rj'K 
puzzle. There was a similar fate for the electroweak sd — > 7*7* — > up contri- 
bution to the short- distance K — > fip amplitude [33]. There is an additional 
experience in the analogous off-shell photon case [14], where the corresponding 
QED quark triangle overshoots the measured value of the related form factor. 

In our case, the smallness of the anomaly contribution to the b — > srj' can be 
ascribed to the depletion of the QCD off-shell triangle vertex. Accordingly, the 
expectation of [8], as we understand it, turns out not to be fullfiled, and we 
seem to be in compliance with the conclusions by [10] that the singlet penguins 
do not do the job and that a combination of several effects within the Standard 
model would be necessary. In addition, our observation of the absence of the 
suppression in the electroweak vertex, when we departure from truly soft glu- 
ons, opens a window for some other contributions which resemble the QCD 
anomaly. For example, ref. [34] advocated that the B — > rj'K puzzle could 
be explained by the additional complicated non-perturbative quark-gluon in- 
teractions related to the anomaly. However, at present these interactions are 
merely parametrized by a phenomenological coupling. Such an extra piece to 
the usual effective Hamiltonian, appearing due to non-perturbative aspects of 
the QCD anomaly, has to be justified yet. We intend to come back to a more 
complete treatment of the B — > rj'K amplitude in a forthcoming paper. 
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gives the ratio 
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